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1. Introduction 

Superfluid 3 He is one of the greatest discoveries in physics [1-3]. It provides a system 
free of impurities and well described by a mean field theory of the Ginsburg-Landau 
type. The development of cooling technologies and detecting devices allows the su- 
perfluid phase not only to be prepared in the laboratory, but its many properties and 
rich phase structure can be carefully measured. This has built a fruitful interaction 
between theory and the experiment. Theoretical predictions can be frequently checked 
by experiments; interesting experimental results also constantly inspire the advances 
in the theory. 
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It is quite intriguing then, that one of the largest disagreements between a the- 
oretical prediction and experimental results has been sitting in this field for over 30 
years. The first-order phase transition between the ABM phase [4, 5] (the A phase) and 
the BW phase [6] (the B phase) has been observed in the laboratory frequently; yet 
the theoretical lifetime of the A phase before the transition (due to standard thermal 
and/or quantum effect) is calculated to be much longer than the age of our universe. 
This is the famous A-B transition puzzle. One may classify the possible solutions into 
3 types : 

(1) Outside influence, such as cosmic rays hitting the sample, as proposed in the baked 
Alaska model [7]; 

(2) Quantum phenomenon, such as resonant tunneling recently proposed [8]; 

(3) Classical transitions. 

Each of these explanations has its own distinctive predictions that can be explicitly 
checked experimentally. 

In this paper, we will not attempt to provide an overview of all proposals. It is after 
all a long standing puzzle in a well-studied system. Our first step is to closely examine 
the status of the baked Alaska model [7,9], which suggests that the transitions are 
triggered by incident radiation. It is the most well-known proposal, and it is regarded 
by some to be the established answer already. At the very least, the experimental 
data [10-14] seem to imply the following almost universally accepted conclusion: The 
baked Alaska model plays a significant role in the A-B transition, but there 
can be competing/parallel effects in certain situations [13]. 

However, we will point out that the same set of data admits a very different in- 
terpretation, one which is at least equally plausible. Both radiation and "compet- 
ing/parallel effects" are only the first step triggering the transitions; Neither 
is the final control mechanism. Of course, this statement is not very interesting 
unless we can propose the final control mechanism, which we will do in this paper. 

Since the 3 He system is the closest known laboratory system that mimics the cosmic 
landscape inspired by string theory [15-17] , any new discovery here may improve 
our understanding of the origin of our universe. It is also possible that some of the 
theoretical advances in the cosmic landscape could be applied here as well [8]. In 
this paper, we will consider the effect of classical transitions [18-21]. It is a recently 
realized phenomenon that first order phase transitions can happen through domain wall 
collisions, even if the colliding domain walls were nowhere near the new phase. This 
new understanding of classical dynamics may be the key to explaining the fast A-B 
phase transitions in superfluid 3 He. 

Based on this recent understanding, we will present the following new proposal. 
Because of the spin and orbital angular momentum of 3 He pairs, typical lab samples of 
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superfluid 3 He have domains and hence intrinsic domain wall networks. There are many 
ways to disturb the network, including incoming radiation and the cooling process. 
However, first-order phase transitions only happen when the disturbances make domain 
walls collide (which can easily happen when the domains move and/or the domain 
walls oscillate). The resulting field (order parameter) excursion following a collision 
allows the quick completion of a first-order phase transition that would have taken 
an exponentially long time otherwise. This means classical transition can be the final 
control mechanism. This proposal yields surprisingly good qualitative agreements with 
many interesting superfluid 3 He observations that cannot be explained by the baked 
Alaska model. 

Furthermore, classical transitions follow a kinematic rule of field excursions, which 
makes predictions very different from models relying on energetics such as the baked 
Alaska. Whether a collision can induce the appropriate field excursion into the B 
phase is controlled by the phase space structure through a single parameter, 8 [22], the 
strong coupling correction from spin fluctuations (paramagnon exchanges) [23]. With 
the appropriate field excursion, the transition will happen if the network disturbance is 
strong enough. On the other hand, without the appropriate field excursion, disturbing 
it more violently will not help at all. This makes a strong falsifiable prediction: points 
on the curve S(P,T) = const, in the P — T diagram should all together be allowed 
or disallowed to transit. Our proposal will be ruled out if transitions can happen on 
(Pi,Tx), but the same or stronger disturbances cannot trigger transitions for (P2,T 2 ) 
on the same curve. 

We will present our findings in the following order. First, in Sec. 2 we give a 
simple review of the mean field theory description of superfluid 3 He, followed by an 
introduction to the A-B transition puzzle. We then will describe the baked Alaska 
model and examine the major positive, negative and neutral evidence for this model. 
We will argue that this evidence should lead one to seek an alternative model, especially 
in light of the theory of classical transitions. 

In Sec. 3 we review the development of classical transitions and briefly explain the 
dynamics. We also made necessary generalizations in order to apply this phenomenon 
to superfluid 3 He — a non-relativistic system of interacting vector fields. 

Finally in Sec. 4, we provide a complete description of our proposal to solve the 
A-B transition puzzle. We calculate the properties relevant for classical transitions in 
this system, i.e., the phase space structure. Details of this calculation are found in 
Appendix A. We list several prior observations which our proposal can explain, and 
provide a definite prediction. 
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2. The Helium 3 A-B Transition 



2.1 Introducing the Puzzle 

The 3 He superfluid phases are BCS states similar to conventional superconductors, 
except that their Cooper pairs have total spin s = 1 and angular momentum 1 = 1. 
This simple difference leads to a much richer structure. 

The superfluid phase is usually described by the order parameter, A ai , a 3 x 3 
complex matrix. Note that it is not a tensor field, since a is a spin index while i 
describes orbital angular momentum. They cannot directly be contracted with each 
other. Spin and angular angular momentum can rotate independently, so the overall 
symmetry group is SO(3) s x SO(3)i x U(l). 

The quadratic kinetic terms are 
i (in 

K= WT>Te ( 9tA ^ A - + cldiKidA*j + c 2 T e ijk dAl k eurndA am ) . (2.1) 

Up to the 4th order, the potential includes all possible contractions: 
\dn ( T \ 

V = 2 de \T C ~ 1 ) A * mAm + ^ lA ^ A ^ A « A « + ^AVA^A^-A^- 

+ ftAVAJAtfAfc- + AAVA^A^ + fcKAeApiKj ■ (2-2) 
Here dn/de is the density of states, T c is the cricitcal temperature, and /%'s are 
parameters we can calcuculate from the theory. Since a is the spin index, only i can 
contract with spacetime derivatives. The orbital angular momentum vector in A Q j 
is just like a vector field, and has transverse and longitudinal velocities ct and cl, 
respectively. 

Above T c , the system is in the normal Fermi liquid state with (A ai ) = 0. For 
T < T c , the symmetry is spontaneously broken, similar to behavior in a standard 
(-m 2 2 + A0 4 ) potential. Here the rich symmetry of 5*0(3),, x SO(3)i x U(l) allows 
two different types of condensation. One of them is the A-phase [4,5], usually described 
by the matrix 

/K0\ 

A(A) — Aa I J . (2.3) 
\0 0/ 

We can see that an 50(3)/ rotation along the z axis is degenerate with a complex 
phase. 

The other condensate type is the B-phase [6], usually denoted 

/100\ 

A(B) — Ab I 1 J . (2.4) 

loo 1/ 
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A simultaneous rotation on both 50(3) 's leaves this "isotropic" phase invariant. 

During the cooling process, the instability into the A-phase actually occurs first. 
Its free energy also remains lower than that of the B-phase down to some Tab < T c . 
There is no way to execute a "super fast" cooling to bypass the above temperature 
range where the A-phase is dynamically preferred. Therefore, the superfluid phase 
transition always goes to the A-phase first in all laboratory experiments. 

Below Tab, the A-phase becomes metastable, the B-phase becomes the true vac- 
uum, and they are separated by a free energy barrier. The superfluid may then undergo 
a first-order phase transition — nucleating a critical bubble of the B-phase that grows to 
convert the entire sample. The tension a and the free energy difference SV between two 
phases are not only theoretically calculated [22,24], but also experimentally measured 
to a good accuracy [25]. We can use these parameters to calculate the critical radius 
R c and the tunneling rate T. 1 

a 

R c = (2 or 3)— ~ 1/im , 
T ~ exp [-S instanton ] < 10- 20000 . (2.5) 

Although a critical bubble can be easily contained in a lab sample, the lifetime before 
this decay, regardless of the dimensionful prefactor in T, is longer than the age of the 
universe. 

Strangely, the A-B phase transition almost always is found to happen when the 
sample is supercooled significantly below Tab- It happens in a few minutes or a few 
hours. In a different system, this kind of discrepancy might be explained by impurities 
or boundary effects. However 3 He superfluid has no impurities. Only 4 He can dissolve it 
in and can be controlled to less than one molecule per mole in the sample. The boundary 
effect actually stabilizes the A-phase [27], as the boundary defines a direction, while 
the B-phase is isotropic 2 . 

This is the famous A-B transition puzzle in superfluid 3 He. 

2.2 baked Alaska Model 

One exotic, but seemingly viable explanation after eliminations, is the baked Alaska 
model [7]. It was proposed that the transition is facilitated by cosmic rays. As highly 
boosted muons pass through the sample, secondary electrons are produced. Those that 
stop do so very abruptly, and thereby act as point sources to locally heat up the system. 

1 The (2 or 3) is the difference between a thermal instanton and a quantum instanton. Their 
rates are both unobservably small and we are quoting the relatively higher one, i.e., the one for zero 
temperature [26]. 

2 Recently, it was argued that the boundary effect, although favoring the A-phase over the B-phase, 
cannot be totally excluded from the discussion [28]. 
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Note that this cannot be thought of as a conventional thermal activation. Ther- 
modynamically, a hot spot cools down by diffusion and the cooling rate is controlled 
by a thermal scale. As mentioned earlier, this process most likely leads to the A-phase 
again. Additionally, during diffusion the exterior is cooler than the interior and so 
always condenses first to whichever phase it is exposed to, which is A-phase. 

Fortunately, as pointed out in [7], the energy released by an electron goes to quasi- 
particles with mean free path ~ 2/im, which is > R c . Therefore, for physics at the scale 
R c , it is not a conventional thermal system. The best description in this situation is 
(still) unknown, so Leggett made a conjecture. 

The baked Alaska conjecture is the following. Since most of the energy is carried 
by quasiparticles with mean free path around 2/xm, below such length scale, there is no 
diffusion and it is appropriate to picture the free motion of quasiparticles as a hot and 
expanding shell. The purpose of this hot shell is twofold. It heats up and then shields 
the interior from the ambient A-phase region to allow an independent cooling. Also, 
regions swept through by this hot shell effectively go through a rapid cooling, which 
avoids the usual cooling constraint and has a chance to settle into the B-phase. Since 
R c < 2/im, it is possible to have regions in the B-phase large enough to expand even 
after the hot shell diffuses away. 

Due to the lack of viable alternatives, this proposal has remained the center of 
focus for over 30 years. Before moving on to examine the evidence, we would like to 
emphasize some conceptual hurdles faced by this portion of the model. It describes 
regions undergoing "heating" and "cooling" on scales below the thermalization scale. 
Namely, it applies thermal intuition on potentially non-thermal dynamics. There might 
be situations where it is a good description, but some skepticism is warranted. Addi- 
tionally, in comparison with data described in Sec. 2. 3, it is a problematic portion of 
the baked Alaska model. 

As an interesting si denote, in Sec. 3 we will mention that during the development 
of the classical transition theory, a similar conjecture was made and later proved to be 
insufficient to provide the whole answer. Together with the experimental data, it is a 
sign that the baked Alaska model could at least be improved. 

2.3 Examining the Evidence 

Our analysis here closely follows [13] which analyzed old data from A-B phase transition 
experiments. It also summarizes the results from a series of more recent experiments 
at Stanford [11,12] which were designed to test the baked Alaska model. 

There are two major positive pieces of evidence. 
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• PI. The nucleation rate becomes significantly higher if the samples are exposed 
to a radiation source. The difference is reasonably proportional to the strength 
of radiations. 

• P2. The temperature dependence is the same (up to an overall scaling) with or 
without the presence of the radiation source. 

There are also three major negative pieces of evidence. 

• Nl. The theoretical prediction has lnT oc —R™ with 3 < n < 5 [9], but the best 
fit to data is n = 3/2, and is unlikely to be n > 2. 

• N2. The measured temperature dependence stays the same when the dominant 
incoming radiation is switched from electrons to neutrons. Theoretically, there 
should be a difference. 

Finally, there is some evidence pointing toward competing/parallel mechanisms. In 
many occasions people found transitions uncorrelated with incoming radiation. 

• CI. A direct search of cosmic ray effect detected no correlation between a cosmic 
ray event and the A-B transition in a superfluid He-3 sample [29]. 

• C2. Transition rates are significantly higher in containers with boundary rough- 
ness (on the length scale fim). 

• C3. Different cooling rates can increase or decrease the transition rate under 
different circumstances. 

• C4. Mechanical disturbances usually cannot trigger the transition, except for 
samples with Si wafers (introduced to create surface roughness) on the brink of 
a transition [14]. 

• C5. At strong magnetic field and low temperature, the magnetic field seems to 
only affect the rate through its effect on the free energy difference AV. However 
in high T/low H, it behaves differently. 

Indeed (PI) establishes a solid correlation between radiation and the phase transi- 
tions. The easiest interpretation of (P2) is that without adding a radiation source, the 
transitions are triggered by the background radiations through the same baked Alaska 
mechanism. However, even if we only consider positive evidence, it does not exclude 
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the possibility that the background trigger is not due to radiation, as long as there is 
a final control mechanism providing the unique temperature dependence. 

Furthermore, (Nl) and (N2) clearly suggest that the baked Alaska model requires 
some modifications. The quantitative prediction that disagrees with data comes directly 
from the baked Alaska conjecture. Together with the theoretical concern mentioned in 
Sec. 2. 2, it is reasonable to seek out alternatives, perhaps in the form of this final control 
mechanism. In fact, according to [13], "While one might argue that radiation could be 
causing nucleation through some other mechanism, it is difficult to imagine another 
and none has been proposed in the 18 years that this problem has been outstanding. ", 
the lack of alternatives seems to be one of the main reasons to keep the baked Alaska 
model in place. 

Thanks to the discovery classical transitions, we do have such an alternative. We 
believe it is worth exploring, especially if it can also explain some of (CI) to (C5). We 
will argue that is indeed the case. 

3. Classical Transitions 

Classical transition is a way to achieve the result of a first-order phase transition through 
purely classical behaviors. Consider a simple Lagrangian as in [19], 



where V is a potential with three local minima, as shown in Fig. 1. It was observed in 
high resolution lattice simulations that when two domain walls between vacua A and 
B collide, the collision region can sometimes go to vacuum C as in Fig. 2. 



Note that although we started by slamming two domain walls together instead of 
being hit by a quantum of radiation, the second half of the induced phase transition is 



c = l(dcf>) 2 -v(4>), 




V(0 




Figure 1: A typical potential allowing classical transitions. 
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Figure 2: Lattice simulation showing a classical transition. The bulk value is false vacuum 
(j)A- Initially there were two bubbles of 4>b- Their collision results in a transition to 4>C- 

very similar to the baked Alaska model. We generate two energetic "shells" to protect 
the middle region, give it a chance to settle into a different minimum, and guarantee 
its subsequence growth. 

Not surprisingly, a very similar conjecture was made to explain this phenomenon. 
Originally in [19] this was described as an "energetic" effect. The highly boosted 
domain walls provide enough energy during the collision to push the field both above 
the potential barrier and toward a field excursion. This predicted that higher energy 
can overcome a higher barrier and reach a 0c further away. 

The prediction was soon realized to be incomplete. In [20] it was shown that a 
threshold energy is needed to cross the barrier, but adding more energy does not push 
the field further. The transition only happens when (<fi c — 0s) < (0b ~ 0a)- 

We think there is a general lesson learned here, which should also apply to the 
baked Alaska conjecture. The collision happens in a very short time scale that our 
intuition of "available energy" does not work, just as the intuition of "cooling" does 
not work in a non-thermal system. A microscopic description is needed, and in this 
case it is the kinematics of the field that controls the transition. 

3.1 The Mechanism 

The classical transition behavior comes from the equation of motion. 



Given two degenerate minima in V and a barrier between them, we can find a 
static interpolating solution f(x) that solves 



with /(— oo) = 4>l and f(oo) = 0#. It looks like a domain wall in the middle where the 
field obviously changes, and stays exponentially close to either vacuum on each side. 
We can arbitrarily boost this solution, /( f=? ) will also be a solution of Eq. (3.2). 



<9 2 <9 2 



dV(<f>) 
d(j) 




dt 2 dx 2 




(3.3) 
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When the vacua are not degenerate, the domain wall will feel a pressure difference 
and a constant acceleration. In the thin wall limit, at every instant the accelerating 
solution is approximated by a boosted solution. 

The behavior of Eq. (3.2) in these highly boosted solutions is crucial to understand- 
ing classical transitions. In the vacuum region it is trivially = 0. Within a highly 
boosted domain wall, both terms on the left hand side are enhanced by a large factor, 
7 2 = 1/(1 — v 2 ). So the zeroth order solution is to ensure they cancel each other, 

d 2 cj) d 2 cj) 

— ~ (3 41 

dt 2 dx 2 ' 1 J 

which is a free field equation. The leftover of this approximate cancellation have to 
agree with — 9 V ^ , but that is akin to next order in perturbation theory. The expansion 
parameter is I/7 2 . 

Now consider three regions in three different vacua, with two highly boosted domain 
walls moving toward each other. Before they collide, the solution looks like 

a j\ 1 r f-X + Vt\ ( X + Vt 

4>{x, t) = 4> M + h n = + f R 



V 



~ <Pri + (0l - <f>M )Q{-x +t) + (<f) R - <f) M )Q{x + t) , (3.5) 

where <f>L, <Pm, <Pr are the vacuum field values of the left, middle, right regions; fz and 
fa are the domain wall solutions from the middle vacuum to left and right vacua, which 
are almost step functions when highly boosted. If all we have to solve is a free field 
equation, Eq. (3.4), then the middle region initially at 

0(O,t<O)=0 Af (3.6) 

just goes to 

0(0, t>0) = <p M + {<j>L - 0M) + (0i? - M =<PL + ^R~ 0M • (3.7) 

This is why adding more energy will not help at all. A fixed "field excursion" is 
determined by microscopic kinematics. If this field excursion is near a vacuum, then 
Eq. (3.2) continues to be solved trivially, and one finds a classical transition into this 
vacuum. This result directly generalizes to multiple fields with standard kinetic terms, 

for which 

0(0, t<0) = M , 

0(0, t>0) = 0l + 0r-0m • (3.9) 



- 10 - 



3.2 Basin of Attraction 



Eq. (3.7) provides a field value right after the collision. For the evolution afterward, we 
can no longer ignore the potential unless it is exactly another local minimum. Usually 
we vaguely say if the field excursion enters the "basin of attraction" of some minimum, 
it will eventually roll there. 

This is somewhat confusing. As in Fig. 3, one might think we just have to go 
over the top of the energy barrier between <pB and 0c 1 , for example to 0o- It is n °t 
this simple. As <p rolls down, there is a competing effect as shown in Fig. 4. Since 
V(4>o) > V((/)b), the middle region will start to give in to the surrounding region in (ps- 
It may recollapse entirely before reaching <pc- 



V 




Figure 3: Showing the naive basin of attraction spanned by the blue arrow, and the absolute 
basin of attraction spanned by the green arrow. 




Figure 4: The field interpolation between </>o and (f>B goes over a barrier, which behaves 
like a quasi-domain-wall. The pressure difference makes them turn around and may collapse 
the middle region before the field rolls from <fio to 4>c- 

On the other hand, this is the only competing effect. So if V((f)o) < V((/)b), then 
it will always roll down to <fic- This allows us to define two regions: the naive basin 
of attraction (just cross the barrier) and the absolute basin of attraction, V < V((/)b)- 
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The real basin attraction which defines whether a classical transition can happen or 
not, must be somewhere between these two. As the temperature below Tab decreases 
towards zero, the size of the actual basin of attraction is expected to grow. 

3.3 Sound Speeds 

In Sec. 3.1 we used the term "highly boosted", which may seem only applicable in 
relativistic systems. That is not true. Adding a sound speed to Eq. (3.2), 

dv{<i>) 

dt 2 s dx 2 " d<J) ' { ' 

is equivalent to rescaling the x coordinate. Alternatively, we can define 7 = (1 — 
v 2 /c 2 ) -1 / 2 , the entire arguement goes through the same way. The domain wall becomes 
highly boosted when it approaches the sound speed. 

The real problem is, in condensed matter systems, we can have multiple sound 
speeds. For example in a system with a vector field, there are the transverse modes 
and the longitudinal mode, with a/3ct = cl- 

}2 ± n2 



d 2 (f) 2 d 2 <p dV((j>,ip) 
dt 2 T dx 2 d(f> 
d 2 i) _ 2 d 2 ip _dV((f),tlj) 
~di 2 ~ CL d^ 2 ~ Ihp 

Now what happens to domain walls in this system? 
First of all, the static domain walls are easy. 

d 2 (c T (j)) dV(<f),ip) 



(3.H) 



dx 2 d(c T (f>) 
d 2 M) dV{M) 



dx 2 d(cLip) ^ 

Simply rescale the sound speeds into the fields definition, = ct<P, if) = cl4>, and it is 
equivalent to a standard multifield problem 3 . 

For a domain wall moving in a constant velocity, we can put in the anzatz (f>(x — vt), 
if)(x — vt) into Eq. (3.11) and get 

( „ 2 _ 4¥ , = _W) (313) 

3 Note that this means the structure (therefore the tension) of a domain wall depends on its ori- 
entation. It is worth further study but quite unlikely to dramatically change the exponentially small 
transition rate. 
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For v 2 < Cy, we can similarly rescale by = \j c T — v 2 </>, ^ = v/c| — u 2 ^, and get to 
a standard multifield problem. But note that the potential seen by the rescaled fields 
depends on the velocity, and therefore the path in the field space (0, ip) depends on the 
velocity, as shown in Fig. 5. For example in the limit v 2 — > c 2 -, the equation of motion 
for (f> is just |^ = 0. This means we should globally minimize V in the <p direction even 
if it leads to a discontinuity. 




(a) (b) 

Figure 5: A domain wall between two vacua at (1,1) and (-1,-1). The left figure is the field 
profile of a moving domain wall. The red curve is the transverse mode </>, which is narrower 
(more boosted) than the longitudinal mode ip (blue curve). The right figure is the domain 
wall trajectory in the field space. The blue path is the same moving domain wall as in the 
left figure. The red path is the stationary domain wall. When v — > ct, the path will become 
a step function. The <j) interpolation becomes very narrow but tp does not. 



How about v 2 > c\l First of all, if the domain wall starts to accelerate from zero 
velocity, the gradient energy will diverge as it approaches or- Even if we try to setup 
an initial condition with something faster than cr, Eq. (3.13) implies an instability in 
the small perturbations of <p. When the domain wall moves faster than ct, it emits 
perturbations traveling at ct- These are no different from Cherenkov radiation! They 
pose a large energy cost, causing the domain wall to quickly slow below ct- As confirmed 
by simulations 4 , this sound speed barrier keeps the velocity of an accelerating domain 
wall at Ct- Only the transverse mode in the domain wall profile becomes highly 
boosted like a step function. The effective boost for the longitudinal mode ip remains 
of order one, 7 = (1 - 4/c£)~ 1/2 = y/3/2. 

Although there are very special potentials (for example Sine-Gordon) where classi- 
cal transitions happen at an arbitrarily low boost, in general 7 = \/S/2 is not enough. 

4 We have only done simple 1 + ID simultions in Mathematica, which is sufficient for our purpose. 
The vector field system itself might worth further study with higher resolution simulations. 
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Here we have a rather peculiar situation — a half classical transition. Indeed classically 
transits but ip does not. 



0(0, t > 0) = L + (f) R - 4> M , 

^(0, t > 0) = p [(lp L - ll) M ) + (ipR - ^m)] + 

= P{^L + i/>r) + (1 - 2p)tfj M . 



(3.14) 
(3.15) 



where p is a factor between and 1 which depends on various details and can fortunately 
be neglected in this paper. 

Note that this analysis is valid for a (1 + 1) dimensional collision, namely a head-on 
collision. In a relativistic system we can always boost to this frame, but in a condensed 
matter system we cannot. We will then assume that there will occur nearly head-on 
collisions in the system for which our analysis is valid, and this can be taken as an 
approximation (or an upper bound) for non-head-on collisions. 5 

4. Proposing a New Answer to the A-B Transition Puzzle 

Having established that classical transitions are part of the field dynamics in this sys- 
tem, we are ready for a new proposal to explain the A-B phase transition puzzle. Since 
we only use the standard equations of motion, our proposal does not involve any novel 
(unconfirmed) dynamics. We will show how the superfluid sample allows domain wall 
collisions to take place frequently. Then we show these collisions can lead to classical 
transitions. 

4.1 A Dynamical Domain Wall Network 

The A-phase 3 He superfluid is a nearly ferromagnetic system. While cooling from the 
normal phase over a time scale r ~ 10 3 — 10 4 s [30], the second-order transition into 
the A-phase has the initial coherence length £ ~ 10~ 2 cm [31]. This is the typical size 
of regions, though all in the A-phase degenerate valley, having their order parameter 
A ai pointing in different directions. These regions are separated from each other by 
domain walls whose thickness is the "healing length" from the dipole alignment effect, 



One may worry that unlike a solid magnet, for which the spins are frozen after 
settling into these regions, the fluid nature can erase such structure. For example, if 
the degeneracy in the A-phase is further broken, the region with the lowest free energy 
might expand and eliminate the others. We think that is unlikely for two reasons. 

5 Actually, part of the results in Fig. 2 are circumstantial evidence showing that non-head-on colli- 
sions lead to the same transition. 



^ ~ 8/rni. [3] 
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The strongest breaking effects come from the container boundary and external 
magnetic fields. Near the boundary, the orbital angular momentum prefers to align 
normally to the surface, and different parts of the boundary will not agree with each 
other. In particular, if the surface roughness has a length scale between £ and it will 
actually pin down the domain wall structure near the boundary. In the middle region 
away from the boundaries, the magnetic moment prefers to align with the external 
magnetic field (if it is strong enough). Note that these are vector effects, so there are 
still degeneracies remaining. There is really no single orientation within the A-phase 
that dominates others through out the sample. 

Second, there is an interesting hierarchy in this domain wall network. The domain 
walls intersect each other in string-like and point-like defects (vortices). Unlike the 
domain walls which are sustained by the small bending energy along the A-phase valley, 
these lower dimensional defects are generally forced to go through A a j = 0, which is a 
peak in the potential. In other words, they are much heavier. In the time scale that the 
domain walls are dynamical, we can treat the strings and vortices as stationary. While 
the domain walls may bend responding to further degeneracy breaking, the strings and 
the vortices might be enough to pin down the overall network structure. 

Therefore, a lab sample is never homogeneously in one particular configuration of 
the A-phase. It consists of many sub-domains of various sub-phases of A, randomly 
formed during the second-order transition from the normal phase. The domain walls 
between them stay in an equilibrium configuration when the sample is not disturbed. 
Once cooling has started, or responding to incoming radiation, the system will be tem- 
porarily out of equilibrium, and domain walls can move to collide with each other 6 . The 
stochastic behavior comes from the repeated collisions between different combinations 
of random sub-phases, as depicted in Fig.6. Each of these collisions has a chance to 
induce a field excursion into the B-phase basin of attraction. 

Superfuild 3 He is a complicated system. There might be more than one dominant 
effect disturbing the domain wall network 7 , and the detailed dynamics of a disturbed 
network are not easy to picture 8 . Fortunately, classical transitions may allow us to 

6 Onc can also imagine nuclcation of sub-phase bubbles, which will expand and eventually collide 
with existing domain walls. The exponentially small tunneling rate in Sec. 2 is only for an A-B transi- 
tion. The transition rate between sub-phases in A has not been estimated and need not be small. We 
prefer to use existing domain walls as our core mechanism for it matches past observations better. 

7 We believe that in experiments with very smooth boundaries [11,12], radiation has the dominant 
effect. In [10], we believe the cooling process and the boundary roughness are the main driving forces 
for domain wall motions. 

8 For example, the strings and vortices, though relatively stationary, can still annihilate with each 
other, and might suddenly "pop" when the surrounding domain walls accumulate a strong stress on 
them. This will subsequently result in further domain wall motion. 
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ignore these details. As long as many collisions happen above the energy threshold, 
the final control mechanism will be the field excursions. 




Figure 6: The double cone is a cartoon picture of the superfluid 3 He phase space structure. 
The middle apex is the symmetric normal phase N. The top and the bottom circles are 
the spontaneous broken phases A and B. Three sub-phases of A (red dots) involved in a 
domain wall collision can lead to a field excursion as the blue arrows demonstrated. This field 
excursion can roll down to the B-phase and complete an A-B transition. 



4.2 Criteria for Classical Transitions 

As mentioned in Sec. 3, a domain wall collision needs to pass a certain energy threshold 
and induce an appropriate field excursion to complete a classical transition. Here we 
will show that both conditions are satisfied in superfluid 3 He. 

The simplest way to think about the energy threshold was discussed in [21]. We can 
just look at Fig. 4 and imagine the collision as a scattering problem of domain walls. Two 
incoming domain walls between sub-phases of A need to become two outgoing domain 
walls between phase A and phase B 9 . gab comes directly from an energy barrier in the 
mean field theory. The same mean field theory tells you that any 2 sub-phases of A are 
connected in a flat direction. Thus, their rest tensions usually satisfy the hierarchy 



This is particularly useful when the incoming domain walls are lighter than the outgoing ones, 
since it guarantees a large boost. For the opposite case it will be more subtle. Fortunately our system 
is the easier case. 
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Therefore, the incoming domain walls need a large kinetic energy to make up the 
difference. 

The kinetic energy comes from the potential difference between sub-phases, Vij = 
\V(Ai) — V(Aj)|. Their precise values are beyond our current knowledge, but since 
boundary effects can totally suppress the A-B transition [27], some V,j is comparable 
to Vab = [VX-A) — V(B)]. In that case, the length scale we need to accelerate the 
domain wall so that its kinetic energy approaches o~ab-, is roughly the critical bubble 
size R c ~ 1/iin. The typical size of a domain in one sub-phase of A is £ ~ lCT 2 cm 
as mentioned earlier. This means we do have enough kinetic energy down to Vy ~ 
Vab /100. It is quite plausible for many collisions to pass the energy threshold. 

As for the field excursions, we put the detailed analysis in Appendix A and only 
summarize the relevant results here. Sub-phase collisions in A never take the order 
parameter A a j directly to the B-phase. This is a desired property. Unlike a nucleated 
bubble, a B-phase region created by a classical transition will be nearly planar and 
surrounded by domain walls still moving apart. Most likely it will prevail as long as 
T < Tab- But experimentally the transition is known to only occur at an even lower 
temperature threshold T ~ 0.7Tab- So if a collision directly created a B region, our 
proposal would have already been ruled out. 

Using the behavior of the naive and the absolution basins of attraction, Appendix 
A also shows that the field excursion starts to go into the B-phase basin of attraction 
at some finite temperature, and moves in deeper at lower temperatures. This means 
classical transitions is a possible candidate to solve the A-B transition puzzle. 

4.3 Agreement with Past Experiments 

Qualitatively, this new proposal of domain wall network classical transitions can agree 
with numerous strange observations made in past experiments, even those ones seem- 
ingly contradicting each other. Many of them, if we stick with the baked Alaska model, 
would have to be explained by competing/parallel mechanisms. 

• Transitions are induced by radiations, but do not agree with the explicit depen- 
dence given by the baked Alaska model. Because incoming radiation only disturb 
the equilibrium and allows further domain wall collisions. 

• Transitions start at T ~ 0.7Tab instead of Tab- This could be the temperature 
for which the possible field excursions start to enter the B basin of attraction. 

• Transitions happen only during cooling but not heating. While bringing a super- 
cooled sample back to a higher temperature, the induced domain wall motion is 
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almost a reverse process of that during cooling. If the involved sub-phases colli- 
sions did not lead to the B-phase basin of attraction during cooling, neither will 
they during heating. This explains why temporarily heating a tube segment in 
the Stanford experiments failed to induce a transition. This also explains why one 
particular sample in [10] did not go to the B-phase at all. Because the random 
A sub-phases in its domain wall network failed to contain a correct combination. 
Within the temperature range they tried, none of the domain wall collisions can 
induce an appropriate field excursion to go to the B-phase. 

• The A-phase is more stable with a smoother container boundary. Because a 
rougher boundary encourages more sub-domains, and may result in nonuniform 
cooling. Both effects increase the chance of collisions between different sub- 
phases. 

• Transitions are generally not induced by mechanical disturbances nor high reso- 
nant RF. Bulk disturbances shift the entire system instead of affecting the struc- 
ture of the domain wall network. To possibly do that, the disturbance needs 
a wavelength comparable to the typical size of a sub-domain ~ lCT 2 cm. The 
frequency needs to be ~ 5 x 10 5 Hz for sound waves and ~ 3 x 10 12 Hz for EM 
waves. 

• As an exception to the previous observation, samples with Si wafers on the brink 
of transitions does immediately go to the B-phase after a gentle tap. [14] Si wafers 
have surface roughness which pins down some of the domain wall network. So 
the relative motion between the Si wafers and the bulk sample does disturb the 
network equilibrium. 

• Fast cooling from T < T c allows the A-phase to reach a lower temperature in 
the Stanford experiments. A shorter process allows less combinations of A sub- 
phases collisions to be tried. Note that fast cooling may allow the domain walls to 
accelerate more as it disturbs the equilibrium faster. However, as long as most of 
the collisions happen above the threshold energy, exploring field space excursions 
is the most important factor in the transition rate. 

• Contrary to the previous point, the Los Alamos group (Boyd and Swift 1993, 
unpublished but results cited in [13]) found that fast cooling from T > T c makes 
the transition easier. The cruicial difference is here the fast cooling started above 
T c in the normal phase. This leads to smaller and more sub-domains forming 
during the second-order N-A transition, which allows more combinations of sub- 
phases collisions to be tried. 
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• Again contrary to the previous two points, Hakonen et al. [10] found that the 
cooling rate has no effect on the transition rate. The container surface in [10] is 
much rougher than the previous two experiments. It leads to many sub-domains 
and collisions between them that the transitions happen rapidly at a catastrophic 
line. It is a sign that a transition is guaranteed once the possible field excursion 
enters the B-phase basin of attraction, therefore independent of the cooling rate. 

• Magnetic fields affect the transitions in a complicated way not limited to its effect 
on AV. In addition to changing the potential which affects the basin of attraction, 
magnetic fields also align the sub-phases of A and changes the possible field 
excursions to be explored. 

Although these agreements are encouraging, they are only qualitative and some 
features can also be explained by the resonant tunneling phenomenon. Next we will go 
on to make an explicit prediction that could be tested by future experiments. 

4.4 Prediction and Discussion 

The two important quantities in classical transitions are the energy threshold and the 
field excursions. In reality, their temperature/pressure dependences might complicate 
analysis and stop us from giving a clean prediction. However since the former is a 
threshold, it should not be hard to setup situations where the domain wall network is 
always disturbed strongly enough to allow energetic collisions. 

In this ideal situation, the transition behavior is determined only by the field ex- 
cursions, and therefore the phase space structure. As argued in Appendix A, it only 
depends on external variables through the quantity 5. Therefore, we can claim the 
following falsifiable prediction for our proposal: 

• Varying only temperature and pressure, a 5(P, T) = const curve determins whether 
transitions are allowed or not. 

Of course, this means leaving everything else the same, including identical cooling 
processes through the N — A transition to provide statistically similar domain wall 
networks. We shall leave the details of how to setup such ideal situations to experi- 
mentalists. 

There is not much room for confusion as our prediction is quite unique. It is only 
shared by mechanisms which only depend on the phase space structure. As far as we 
know, the only existing proposal that might have this property is related to Q-balls [32], 
but we are not aware of a specific predictions made in that direction. 

Finally, we should point out a logical by-product of our overall argument in this 
paper. It is possible for the domain wall network and classical transitions to enter the 
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picture in a less dramatic way. First, they provide the "competing/parallel" mechanism 
since our list in Sec. 4. 3 explains all of the (CI) to (C5) in Sec. 2. 3. Next, we have 
learned two possible ways to improve the only weakness in the baked Alaska model — 
the quantitative prediction. 

The first possibility is to keep the expanding shell conjecture as temporarily pro- 
tecting the middle region, but instead of a fictitious "fast cooling", the middle region 
should be described by some kinematic rule of field excursions. So instead of an "ener- 
getic" prediction which does not agree with experiments, the baked Alaska model can 
have a kinematic prediction similar to classical transitions. 

The second possibility is to trigger domain wall collisions with incoming radiation 
but assume the field excursion is mostly satisfied, so the energy threshold becomes 
the control mechanism. Although this is again an "energetic" prediction, disturbing a 
2D domain wall is geometrically different from exciting a hot shell. This can explain 
the lack of difference between electron and neutron radiations, and provide a more 
appropriate value of n in the relation T oc —R™ to fit the data better. 

Both of the above proposals require further study, and their predictions are unlikely 
to be as elegant as our main proposal. Thus, at this stage, we think the easiest way to 
determine the significance of classical transitions is still the prediction we mentioned 
earlier: Transitions are controlled by the phase space structure, therefore S(P,T). 
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A. The Phase Space Structure 
A.l Direct Transitions 

Consider a general SO(3) s x 5*0(3); x U(l) rotation acting on the standard A-phase 
matrix, Eq. (2.3). The general form of all degenerate A sub-phases can be written as 



Here s is a real unit vector for the spin, whereas u and v are real unit external vectors 
orthogonal to each other, u ■ v — 0. There are 5 degrees of freedom here. An overall 
complex phase is not necessary, since it is degenerate with the rotation along u x v. 

First, we would like to find three sub-phases in A to call the left, middle, and right 
regions, and then apply Eq. (3.16) to see if the classical transition goes to the B-phase 
directly. For a head-on collision, we can rotate the domain wall motion into the x-axis. 
Then s a u x and s a v x play the role of the longitudinal mode ip, while s a u y , s a v y , s a u z 
and s a v z are the transverse modes, like <fi. Hence we solve the following set of equations: 



A(A) = A A s(u + iv) . 



(A.l) 



p(ulxSl + u Rx s R ) + (1 - 2p)u Mx SM 



p(vlxS L + v Rx s R ) + (1 - 2p)v M xS M 



U Ly SL + U Ry S R - U M ySM 



ULzSL + Ur z Sr - UMzSM 



VLySL + VRySR - V My SM 



VLzSL + Vr z Sr - V Mz Sm 
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sin 5(1,0,0) , 



sin 9(0, 1,0) , 



sin 9(0, 0,1) , 



cos0(l,O,O) , 



cos 9(0, 1,0) , 



cos 9(0, 0,1) . 



(A.2) 



Note that the left-hand sides are vector equations of spins, so we can still use 
SO(3) s to rotate the right-hand side into the diagonal form of the B-phase. However the 
complex phase 9 cannot be rotated away since the choice of collision axis already breaks 



S0(3)i. Counting the degrees of freedom, we have (5 x 3 + 1) = 16 free parameters, 
but 18 equations, so the system is overconstrained. 

To explicitly see that a direct transition fails, we introduce the vector 

w = v cos 9 — u sin 9 . (A. 3) 

Since v ■ u = and v, u ^ 0, w must be a nonzero vector. This allows us to pair up the 
six equations (A. 2) for the following three: 

p{w L xS L + w Rx s R ) + (1 - 2p)w Mx s M = (0, 0, 0) , (A.4) 
W Ly S L + W Ry S R - W My SM = (0, 0, 0) , (A. 5) 

w Lz s L + w Rz s R - wmzSm = (0, 0, 0) . (A. 6) 

Because the w's are nonzero vectors, it is impossible to make all three equations 
trivial. At least one of them must tell us that sm is just a linear combination of sl 
and sr. This implies that the left-hand sides of the last three equations in Eq. (A. 2) 
are limited to a 2D surface, while the right-hand sides are three mutually orthogonal 
vectors. That is impossible. 

Thus, we can conclude that classical transition cannot bring the field value directly 
to the bottom of the B-phase. The next step we need to consider is the effect from its 
basin of attraction. 

A. 2 Basin of Attraction 

We need to take a closer look at the potential, Eq. (2.2) to understand the basin 
of attraction. Following [22], the coefficients Pi come from the weak coupling BCS 
calculation 

_2/3 1 BCS = /3 2 BCS = /3 3 BCS = /3 4 BCS = -/3 5 BCS , (A.7) 
while the strong coupling correction provides 

/3 2 = /3 2 BCS (l + ^), 
/3s = /3 3 BCS (1 " ^S) , 

/3 5 = /3 5 BCS (l + ^). (A.8) 
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Here S = 8ab — 20/43 is the critical value when Va = Vg . Smaller S makes the B-phase 
more stable. 

Note that the difference between A and B phases A are only in their 4th order 
contraction terms. Thus the ratio between the 2nd order and the 4th order coefficients 
can always be scaled away by redefining the lengh |A a j|. This will not affect the relative 
basin of attractions between A and B. So the only relevant change is through 6, as it 
affects all the 4th order coefficients differently. We can hold everything else fixed and 
only vary 5 to capture all physical effects relevant to the basin of attraction. 

As argued in Sec. 3. 2, the real basin of attraction is hard to delineate, but it must 
lie between the naive and the absolute basins of attraction. 

For the naive basin of attraction, we perform the following numerical search. We 
randomly generate three sub-phases of A to be the left, middle and right regions (5x3 
random variables) and combine them with the classical transition rule, Eq. (3.16). 
Starting from this value, we follow the gradient flow (— VV) (an 18 dimensional vector) 
to see if it rolls down to the B-phase. 

Since < p < 1 is unknown, we tried three values, 0, 1/2, 1, and the results are 
similar. As shown in Fig. 7, decreasing 5 continuously increases the chance for it to 
enter the B-phase. This implies the possible value of A ai does indeed go deeper and 
deeper into the B-phase naive basin of attraction. 10 

As for the absolute basin of attraction, our approach is to first construct the 15- 
dimensional compact subspace (in the 18-dimensional full phase space of A ai ) that is 
reachable by a classical transition. Then we locate the global minimum in this subspace. 
If this restricted global minimum satisfies V g < Va, it must be in the absolute basin of 
attraction of phase B, since it is the only other minimum to roll to. 

Again we tried three values of p and they all have the same behavior. For 5 > 0, 
this compact subspace never intersects the B-phase absolute basin of attraction. But 
zero is a critical point: Once 5 < 0, this compact subspace enters the absolute basin 
of attraction. Note that negative 5 is not a physical choice of parameter, and in fact 
there is a second-order phase transition here whereby the A-phase becomes a saddle 
point. However we are just looking for a statement of continuity. Since the potential 
V and the structure of the B-phase are continuous functions of 5 even through 5 = 0, 
we can conclude that the B-phase absolute basin of attraction expands as 5 decreases 
and eventually touches this 15 dimensional subspace. 

10 It is not a relative effect that somehow the distance between the B-phase and the 15 parameter 
subspace of possible collision results are getting closer. Rather, this is a global effect that the B- 
phase (absolute) basin of attraction is getting larger. We determined this by randomly choosing 10000 
arbitrary A a i instead and tracking to where they roll. The result is just Fig. 7 multiplied by a constant. 
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Figure 7: The number of times A Q j is within the B-phase naive basin of attraction, as a 
function of 5. For each 5 there are 10000 random starting points as possible values reached by 
collisions. 5 = 20/43 corresponds to Tab i n this calculation. The "partial classical transition" 
parameter p = 1/2 in the data set here. For p = and p = 1, the values are roughly lower 
by 25% and the trend remains the same. 



We have established that the naive basin of attraction covers more and more values 
of A a i reached by collisions, and that the absolute basin of attraction approaches and 
eventually touches these values. Therefore, the real basin of attraction, which is always 
between these two, must start to cover these values at some finite S, and likely covers 
more values at smaller S. 
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